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Abstract 

Recoupling coefficients (3nj symbols) are unitary transformations between binary 
coupled eigenstates of N = (n + 1) mutually commuting SU(2) angular momen- 
tum operators. They have been used in a variety of applications in spectroscopy, 
quantum chemistry and nuclear physics and quite recently also in quantum gravity 
and quantum computing. These coefficients, naturally associated to cubic Yutsis 
graphs, share a number of intriguing combinatorial, algebraic, and analytical fea- 
tures that make them fashinating objects to be studied on their own. In this paper 
we develop a bottom-up, systematic procedure for the generation of 3nj from 
3(n — l)j diagrams by resorting to diagrammatical and algebraic methods. We 
provide also a novel approach to the problem of classifying various regimes of 
semiclassical expansions of 3nj coefficients (asymptotic disentangling of 3nj di- 
agrams) for n > 3 by means of combinatorial, analytical and numerical tools. 



Keywords: Quantum theory of angular momentum; Racah-Wigner algebra; 3nj 
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Introduction 



In the quantum theory of angular momentum -mathematically encoded in SU(2) 
representation theory- the study of (re)coupling of pure eigenstates of several, 
mutually commuting angular momentum operators J 1; J 2 , . . . Jtv is considered an 
advanced topic, and as such it is not widely known. As pointed out in the classic 
reference [1J (topic 12) 

We shall deal with what is appropriately called the theory of binary 
coupling of angular momenta, since it is a theory in which one cou- 
ples angular momenta sequentially, in pairs. (...) Accordingly, the 
theory deals with the relationships between different coupling schemes 
-that is between distinct sequences of pairwise couplings. This is the 
content of recoupling theory. 

The subject is a difficult one and the literature is extensive. Several 
monographs {e.g. B2l[3l|4l|5][ll) deal with techniques for implementing 
"graphical" methods (...). 

This paper is much in the spirit of the program stated soon after that, namely 

Our goal is to relate coupling methods to standard results from graph 
theory. The subject divides quite naturally into three parts: 

(A) (the classification of coupling schemes) 

(B) the elementary operations underlying the structure of transforma- 
tion (recoupling) coefficients, or 3nj symbols 

(C) the classification of recoupling coefficients and its relationship to 
the theory of cubic graphs. 

We are going to deal first with (B) and (C) above by formalizing an effective, 
step-by-step method to build up 3nj coefficients of type I, II from 3(n — l)js of 
the same type (note that, by convention, the integer n is related to the N angular 
momenta introduced above by (n + 1) = iV). 

Our second goal is to take advantage of such a recursive construction in the 
search for a novel unifying scheme for addressing semiclassical limits (asymptotic 
expansions) of the 3nj themselves. This task calls into play different techniques, 
both analytical and numerical, and is very important in view of practical applica- 
tions (that have been discussed extensively elsewhere). 

'The last two monographs were published seven years after Biedenharn-Louck text; in par- 
ticular the last one is now the most exhaustive and reliable collection of definitions of angular 
momentum functions and related algebraic and analytic formulas. 
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1 Recursion method for 3nj coefficients: 
general description 

The diagram of any SU(2) 3nj coefficient is a cubic (regular trivalent) graph on 
3n lines and 2n nodes, named Yutsis graph from the first author of 0. The three 
edges stemming from each node are associated to a triad of angular momentum 
quantum numbers satisfying triangle inequalities, i.e. to a Wigner 3j coefficient. 
This property of diagrams reflects the fact that any 3nj can be always written as 
a sum over all magnetic quantum numbers of the product of 2n (suitably chosen) 
3j coefficients (see e.g. O, Ch.9). 

The recursion method for the generation of (all) 3nj diagrams with increasing n 
may be summarized as in [2] p. 65 

(Recursion rule) The diagrams of the 3nj coefficients are obtained 
from the diagrams of the 3(n — l)j by inserting two additional nodes 
on any two lines of the diagram of the 3(n — l)j coefficients and 
joining them together. 

There are some caveats underlying this combinatorial construction: it must gen- 
erate diagrams that are not separable on fewer than four lines and this is achieved 
by requiring that the initial 3nj diagram is separable on no fewer than three lines 
(examples are provided below); moreover, isomorphic configurations should be 
recognized and ruled out||. 

Algebraic expressions ("recurrence formulas") that encode the above combi- 
natorial prescription on diagrams are provided for the 3nj coefficients of the first 
(I) and second (II) type in ((24.8) p.81 and (A.6.14) p. 143), respectively and 



where $ = j x - &i - j n _ x + fc^-i, and 

2 As an aside remark recall that the isomorphism problem in graph theory is a typical NP- 
complete problem, actually a #P problem, even for regular graphs: this means that, given two 
instances of graphs, one can decide if they are isomorphic or not in polynomial time but the general 
solution -the enumeration of all graphs isomorphic to a given one- requires an exponential time 
as the complexity of the graphs grows, i.e. as the number of vertices increases (see e.g. ||6]). 
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where ^ = j\ — k\ +j n -i — k n -i- Here the angular momenta variables {ki, . . . , h, 
. . . , ji, . . . } run over either {0, 1, 2, . . . } or {|, |, . . . } and satisfy suitable trian- 
gle inequalities. The summation variable x is also constrained by triangle rela- 
tions. 



The roles played by the expressions (QQ) and © are multiple and interconnected 

i) at facevalue they provide a straightforward decomposition of a 3nj coeffi- 
cient into a single sum involving one 3(n — l)j of the same type and exacly 
two Wigner 6j symbols. 

(Note however that these decompositions are not the basic ones for compu- 
tational purposes, since it is well known that both types of coefficients may 
be expressed as single sums of products of n 6j symbols [5], Eqns. 1 and 2 
p. 361.) 

ii) from the structural point of view they encode also generalizations of the 
basic Biedenharn-Elliott identity (see the case of 9j generation in section 2 
where these two aspects of (OQ) and © are worked out explicitly). 

iii) taking advantage of symmetries of (any kind of) 2>nj coefficients one can 
equate two different decompositions related by a symmetry of a same 3nj 
to get relations for the associated 3(n — l)j coefficients. For instance the 
relations derived from (OQ) have the following general structure (dropping 
phases and omitting all the arguments not relevant in the present discussion) 



E 



E 





x 



(3) 



A 



Specifying A = |, 1, |, 2, etc. and using explicit forms of the resulting 6j 
symbols, various recursion relations may be obtained (note that the sums 
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on both sides of © reduce to a same number of terms, e.g. for A = | one 
would get x, y — 0, 1). The case of the 9j symbol and associated 5-term 
recursion relation will be worked out in section 2.2. 

The case of the 6j symbol is quite subtle since the well known 3-term recur- 
sion relation in one variable - used in EUl, EL topic 9) as a second order 
difference equation to be analyzed within the WKB framework- actually 
derives directly from © which is nothing but the standard Biedenharn- 
Elliott identity (see also © in the next section). This remark has to do with 
a deep property of hypergeometric polynomials in the Askey scheme BUl: 
the defining difference or differential equation and the recursion relation 
are usually "dual" to each other but the 6j happens to be "self-dual" so that 
the two viewpoints can be used equivalently. 

iv) Decompositions like those given in © and ©, as well as associated recur- 
rence relations, can be also employed in the opposite direction, namely to 
show how a 3nj symbol falls back into a 3(n — l)j when one of its entry is 
set to zero. This well known feature O will be used in remark (i) at the end 
of section 2.2. 

The bottom-up recursion method, fully developed in this paper by resorting to 
both diagrammatical and algebraic tools ©|3l|5l, represents not only a systematic 
procedure for the generation of 3nj's from 3(n — l)j diagrams (section 3), but 
will provide also a novel approach to the classification of the (various regimes 
of) asymptotic disentangling of the coefficients themselves (section 4), improving 
and extending previous results for the 9j found in IfTOl . 

To make the reader familiar with diagrammatic and algebraic tools, in the next 
section we are going to address the case study given by the 9j (1,11). 

2 Genesis of 9 j coefficients and associated diagrams 

The 9j(I) coefficient (the "standard" 9j) can be written as a single sum of three 
6j symbols (notations as in [0, Ch. 10 and App. 12). Note that the following ex- 
pression can be looked at as the first instance of the recursion formula ©, namely 
this coefficient arises as a 3nj(I) for n = 3 starting from the unique 3(n — l)j 
with n = 2, i.e. from the 6j symbol itself 
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where max{|a — b\, \f — e\, \c — d\} < x < min{a + b, f + e, c + d}. 
The associated Yutsis diagram is 



a 




b 



from which the triad structure as well as the symmetries of the coefficient are 
easily inferred (the six triads are associated with 3-valent nodes and correspond 
to columns and rows of the array; symmetries, up to phases, are implemented by 
odd or even permutations of columns or rows and the value of the coefficient does 
not change under transposition). 

On applying type II recurrence formula © one should get a 9j(II) coefficient 
but it can be easily shown that the resulting configuration is a "trivial" one, being 
the product (without sum) of two 6j's sharing a common triad 




The content of this formula is however highly non-trivial since it can be recog- 
nized as the Biedenharn-Elliott identity relating five 6j coefficients. (Recall that 
it represents, together with the orthogonality conditions, the "defining relation" of 
the hypergeometric polynomial of type 4F3 associated with the 6j, see [5] p. 295 
and ®.) 

The associated Yutsis diagram is clearly "separable" on the three lines p, q, r 

P 




r 
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where on the right there appear two Yutsis diagrams of the 6j symbol, each to be 
thought of as a complete quadrilateral whose trivalent nodes correspond to triads 
of angular momentum variables according to the convention 




The two inequivalent ways of applying the Yutsis rule stated at the beginning 
of section 1, namely "inserting two nodes and joining them", generate in this 
simplest case two graphs isomorphic to the 9j(I) (left) and 9j(H) (right) diagrams 




The explicit relationship between such a combinatorial procedure and the content 
of the algebraic recursion formulas © and © is postposed to section 3, where 
the insertion operations will be addressed in the 3nj(l,ll) (any n > 3) cases and 
related to the general expressions given in (OQ) and ©. 

2.1 Combinatorial properties of 6j and 9j graphs 

In this section we are going to discuss properties of the simplest Yutsis graphs 
which will be recovered in an improved way when addressing recursively 3nj 
(1,11) diagrams in section 3. 

a) The diagrams considered so far possess Hamiltonian circuits, as shown be- 
low. 
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Recall that a Hamiltonian path in an undirected graph is a path which visits 
each vertex exactly once. A Hamiltonian cycle, or circuit, is a Hamiltonian path 
that returns to the starting vertex. Determining whether such paths and cycles ex- 
ist in graphs is a typical instance of NP-complete problem [[6]|. 
As will be shown in section 3, both 3nj graphs of type I and II admit by construc- 
tion a Hamiltonian circuit of length 2n for all n, bounding a polygonal region 
("plaquette") with the same number of sides (this is not necessarily true for other 
types of diagrams encountered for n > 4). 

b) Apart from Hamiltonian circuits, the diagrams possess other cycles bound- 
ing characteristic polygonal plaquettes. The drawings below display, for 
each diagram, a typical cycle of this kind (of course, owing to symmetries of 
the coefficients, other cycles of the same shape might have been depicted). 




The number of edges N^A of the "largest", non-Hamiltonian circuit in a 3nj 
diagram depends only on n (there are only triangular plaquettes for n = 2 while 
there appear quadrilaterals for n — 3, pentagons for n = 4 etc. , so that = n+ 
1). The 9j (II) graph actually contains also triangular circuits (not highlighted in 
the above picture) and this feature is related to the "separability" of this particular 
diagram as discussed above. Since the girth of a graph is defined as the length of 
the shortest cycle contained in the graph, the 9j(II) diagram is characterized by 
girth 3 while the 9j(I) has girth 4 = N$L. 

Generally speaking, the presence of cycles with a number of edges strictly less 

(n) 

than iVmax is related to the separability of 3nj coefficients into a product (no sum) 
of lower order coefficients, namely to "trivial" configurations. 
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c) The existence of circuits (Hamiltonian or not) can be assumed as a "guiding 
principle" in the search for asymptotic expansions of the 3nj coefficients 
by letting the edges belonging to the circuit to become "large" while keep- 
ing "quantum" the others. This is what has been already done for the 6j in 
ifTTTl . for the 9j in ifTOll and general results on the so-called phenomenon of 
"disentangling" of networks will be discussed in section 4 below. 

d) By resorting to notions from topological graph theory (see section 3) it is 
possible to classify Yutsis graphs associated with 3nj(L, II) diagrams ac- 
cording to the possibility of "embedding" them on closed surfaces. 

Consider in particular the graphs of the 6j (to be classified as type II!) and of the 
9j (II) (the separable one). They are both embeddable on a 2-sphere since they 
represent (the surfaces of) a tetrahedron and a triangular prism, respectively 



On the other hand, the 9j(I) diagram is a genus-1 or "toroidal" graph, as can be 
inferred by noting that this graph is isomorphic to (the bipartite graph on 6 
nodes, or "utility graph"). To make this point explicit, recall that the torus T 2 can 
be thought of as a rectangle in the plane with opposite edges pairwise identified 
without twisting. The existence of a Hamiltonian circuit connecting all of the six 
nodes is crucial here: the circuit with its nodes is drawn as an horizontal straight 
line inside the rectangle and represents a longitudinal path on the torus (which is 
topologically a circle owing to identification of opposite points on vertical bound- 
aries). The other connections are established in such a way that the connectivity of 
the original diagram is recovered (points to be identified pairwise have the same 
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horizontal coordinates on the upper and lower boundaries). 




In section 3.2 the embedding of both 9j(I) and any other 3nj(l) diagram into 
the projective space will be addressed. 



2.2 Recursion relation for the 9j symbol 

According to remark iii) of section 1, the starting point is the generating formula 
for a type I 3nj coefficient given in (QQ)- It can be used to get an expression 
involving (sums of) products of two 3(n — l)j and four 6j as in ©. The explicit 
form of the latter in the present case reads (0 p. 345; the same in [0 p. 175) 

fa b x 
d e f 



9 h j 



f+f'+9+-J2{-) 2y (2y + l 





Specifying A = 1, a' = a and f = f and using explicit expressions for the 
6j symbols, the above formula gives a 5-term recursion relation for the 9j ([5], 
(10.5.3) p. 347; the same in Q p. 177), namely 

A c+l (ab,fj) fc+i^l A c (ab,fj) f c-i a b ^ f \f c f a , b \ 

-{ f d e \ -\ — i f d e } -f C(a,0,C, f, 7 < f d e > 

(C + 1)(2C + 1) I j 9 hi C (2C+1) I j 9 hi V ' ' ,J,JJ \jgh) 

(7) 

Ad+i(ef, ag) ( e a 6 1 A d (ef, ag) r c a n j c a f> ^ 

where symmetries of the 9j have been used. The functions in front of the 9j 
coefficients are given explicitly by 

A q (pr, st) = [(—p + r + q)(p - r + q)(p + r - q + l)(p + r + q + 1) (8) 

l 

x (-s + t + q)(s - t + q)(s + t - q + \){s + t + g + l)]a 
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C(a,b, c, f,j) 




V{a,d,g,e, f) 




The meaning of these quantities is geometrical, as happens in the 3-term re- 
currence relation for the Qj symbol Q, flTJ (topic 9), 0. (Note however that the 
simplest derivation can be found in |[T2ll . where it was shown that only 6j's where 
one of the entries is 1 are involved, as happens also here for the 9j.) In particu- 
lar, A q (pr, st) represents the squared area of the quadrilateral bounded by p, r, s, t 
and made of two triangles meeting at q, while C and V are cosines of the dihedral 
angles at c and d, respectively (see [fT3l . section III for more details). Each line in 
© closely resembles the structure of the recurrence formula of the 6j, the main 
difference being that here there are two running variables, (c, d) instead of one 
(the other spins play the role of parameters). Notice that the entry labeled by h 
is not active, but we might write down a similar expression involving either the 
pair (c, h) or (d, h), and finally collect all of them in a 7-term relation in the three 
running variables c, d, h, as we are going to sketch below. 

By introducing the shorthand notation 



where % h represents the function to be associated with the third running variable, 
the complete recursion equation can be casted in the form 

A c+ iAf(c+ l,d, h) + A c Af(c- l,d, h) +C c Af(c,d, h) 

= A d+x N{c,d+l,h) +A d N{c,d-l, h) +V c AT(c,d,h) (12) 

= A h+ iN(c,d,h + 1) + A h J\f(c,d,h- 1) +U h Af{c,d,h) 

(the symbols A include the proper normalization factors defined as in ©). 

In the (triple) asymptotic limit (c, d, h) ^> 1 — >■ (x, y, z) the "difference" equation 

above becomes a partial differential equation for the function J\f(x, y, z) 

[A 2 X - 2 cos 9 X + 2] M{x, y, z) = [Aj - 2 cos 6 y + 2] Af{x, y, z) 




V(a,d,g,e,f) = V d ; H{b,e,h,g,j) 



(11) 



[A^ - 2 cos 9 Z + 2] Af(x,y,z) 



(13) 
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where A 2 , Ay, A 2 Z denote second (partial) derivatives with respect to the con- 
tinuous variables x, y, z while 9 X , 9 y , 6 Z are angles whose cosines are defined 
in terms of both x, y, z and the other parameters of the original 9j coefficient. 
The derivation of each member of (fT3l) can be carried out as in [8] and we refer 
the reader to [fT3l for details (see in particular their Eq. 22). However, in the 
present case the proper domain of the function J\f(x, y, z) is the Euclidean space 
R 3 and (A 2 ,, A 2 , A 2 ) = A 2 is the 3-dimensional Laplace operator. Moreover, 
if we equate each member of (fT3l) separately to the same function J\f times a 
constant (to be determined), we would get a set of coupled eigenvalue equations, 
and the allowed values for the separation constant would generate the eigenvalue 
spectrum. Two more remarks are in order 

(i) By setting h = in © the conditions b = e and j = g follow necessarily. 
Then, according to remark iii) and iv) of section 1 , we would see how the 
5-term recursion for the 9j simplifies for the 6j into a 5-term recursion and 
thus, as we have just seen, in a two-variable difference equation which, in 
the asymptotic limit, leads to a two-dimensional coupled partial differential 
problem. Though such kind of procedure seems rather trivial because we 
already know that the 6j satisfied a 3-term recurrence for which the "sepa- 
ration constant" is given, this point illustrates the "dual role" of, say, c and 
d, and the interest of studying the 6j symbol as a function of two variables 
at fixed values of the other four entries (see [fiTTl ). 

(ii) The above procedure, when applied to higher 3nj symbols, can be shown 
to lead to exactly n coupled equations in n variables for symbols of the first 
and second kind, and for highly symmetric ones, such as the 15j of the fifth 
type, whose Yutsis representation is the Petersen graph. Less symmetric 
cases might be dealt with according to the degree of their edge-transitivity. 

More generally, such asymptotic techniques based on recurrence formulas - 
developed here for the 9j- represent the starting point for further developments 
concerning Rodriguez-type formulas (the defining differential relations of fami- 
lies of orthogonal polynomials, see O [fT4l for 1 -variable polynomials), recursion 
relations lfT51l and their relationships with the issue of separability of Schrodinger 
equation in many-body quantum systems, see lfT3l and references therein. Apart 
from a short account in lfl4l (Ch. 5, Sec. 4), such kind of analysis in the case of 
2-variable orthogonal polynomials of hypergeometric type -as the 9j happens to 
be- seems to have been addressed only quite recently in the mathematical litera- 
ture [16]. We argue that results we are going to review in section 4.1 on particular 
asymptotic expansions of the 9j and the novel general scheme for asymptotic dis- 
entangling of 3nj discussed in section 4.2 might contribute to shed light also on 
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the generation of hierarchies of multi-variable special functions of the hypergeo- 
metric type. Work is in progress in this direction. 

3 Genesis of 3nj diagrams 

Looking at Yutsis diagrams as particular families of (unoriented) graphs, namely 
collections of nodes (vertices) and links (edges) connecting pairs of nodes, it is 
worth to introduce a few notions and definitions from topological graph theory 

ma. 

A graph is embeddable on a surface £ if its vertices and edges can be arranged 
on it without any crossing. The "genus" of a graph is the lowest genus of any sur- 
face on which the graph is embedded. Recall that closed and orientable surfaces 
E are topologically classified by their Euler number = 2 — 2g, where the 
genus g is the number of "handles" or "holes"; then x — 2 (g = 0, namely no 
handles) gives the 2-sphere S 2 , x = (g = 1) gives the 2-torus T 2 (whose pre- 
sentation has been shown at the end of section 2.1), while the cases g — 3, 4 . . . 
correspond to double, triple, ... torus S s , all characterized by an even negative x- 
A planar graph is one which can be drawn on the Euclidean plane or, equivalently, 
on the sphere S 2 without any crossing, namely it is has genus or x = 2. 

The issue of the classification of Yutsis graphs is a longstanding problem 
Q3 [HI and in particular such a "topological" viewpoint has been addressed in 
lfi~9l . In this section we are going to complete the characterization of all 3nj(I,II) 
diagrams by resorting to the analysis of the combinatorial and topological contents 
of the "insertion operations" in connection with the recurrence formulas (OQ) (type 
I) and © (type II). Moreover, the existence of Hamiltonian cycles (defined below) 
in any of such diagrams will be the starting point of the analysis of asymptotic dis- 
entangling developed in section 4. Finally, a few remarks on the application of our 
recursive procedure to other types of 3nj diagrams are briefly discussed. 

We start from the analysis of type-II graphs since they share simpler features. 

3.1 3nj(II) diagrams 
Proposition 1. 

Any 3nj diagram of type II is planar (or genus-0). 

Recall from section 2.1 that both the 6j (tetrahedron) and the 9j(II) (triangular 
prism) graphs are easily recognized as planar. Comparing the original diagrams up 
to the 18 j given in it turns out that a convenient and unified way of drawing 
all 3nj(II) is through Schlegel diagrams of prisms based on polygons with an 
increasing number of sides. Below the square (12 j), pentagonal (15j), hexagonal 
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(18j) prisms of this sort are depicted, while the triangular prism associated with 
the 9j(I) has been already shown at the end of section 2.1. 




Other non-Hamiltonian cycles of increasing length appear as well, and it turns out 
that the girth number -the length of the shortest cycle- is exactly 4 for all of these 
graphs (recall that the property of having girth 3, namely the existence of triangu- 
lar circuit(s) implies separability of the diagram into lower order diagrams and as 
such has been already ruled out, cfr. the recursion rule stated at the beginning of 
section 1). 

Proposition 2. 

There exists (up to symmetry) one kind of insertion operation, denoted by 3u, 
that generates the 3nj(II) from the 3(n — l)j(II)for any n > 3 according to the 
recurrence formula ([2]). Moreover J n is topology-preserving, namely all graphs 
derived in such a way stay planar. 

The proof can be carried out by resorting inductively to Yutsis graphical method 
which mimics step-by-step algebraic manipulations to be performed on equa- 
tions. The figures below deal with the generation of the 12 j (IT) from the 9j(H) 
and make it manifest the effective equivalence between the combinatorial recur- 
sion method and the algebraic recurrence formula. 
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In the first figure the recursion rule "add two vertices and join them" (pictorially 

© = ©) is applied on any one of the lateral faces of the triangular prism. 

In the second one the content of the algebraic recurrence formula © for 3nj(ll), 

n = 3, is displayed: the summation over x implies the cancellation of both edges 

labeled by x and nodes which lose an incident edge; in the last step edges with the 

same labels are joined again to give the final closed diagram. 

Note that the two diagrams of the 9j(H) on the left-hand side of each picture are 

isomorphic, as well as the two resulting graphs of the 12j(II) on the right-hand 

sides. 




Combinatorial recursion rule (type II) 




Algebraic recurrence formula 9j(II) — > 12j(II) 

The diagrammatic (box-like) and algebraic forms of the insertion operator 
to be used recursively for the generation of any 3nj(II) diagram or associated 
coefficient are summarized in the following table (where the factor (2x + 1) and 
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a phase have been dropped) 



\a by 



„ [a b xl (a' b' xl 
d , C ^ E *\c d e)\c d f) 



3n ^ 



a/ f \6' 



Noticing that such an insertion can be applied by construction without changing 
the intrinsic topology (any prism with basis an n-gon is changed into a prism with 
basis an (n + l)-gon) the proof of Proposition 2 is completed. ■ 

3.2 3nj(T) diagrams 

Looking at samples of 3nj(I) diagrams as depicted in (2l[3]|, it is easily inferred 
that all of them are "cartwheel" configurations, labeled consistently as shown be- 
low (the arrangement of labels complies with the right-hand side of CQ)) 



Recall that the computation of the Euler number x of a connected configura- 
tion of nodes and links drawn in the plane can be performed by evaluating the 
quantity V — E + F (V = number of vertices, E = number of edges, F = 
number of faces). Here a "face" is any region of the plane bounded by some 
polygonal contour with a certain number of vertices, but graphs with multiple 
edge-crossings, as those in the previous picture, must be suitably rearranged in 
order to get the correct counting. To such a number it must be added 1 to take in 
account the "external" region, namely the portion of the plane not included into 
the graph. 

Proposition 3. 

Any 3nj diagram of type I is a graph with Euler number x — 1. 

The proof proceed by induction on n and will be completed after Proposition 
4. According to the definition given in the introduction to this section, it is worth 
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noting preliminary that an Euler number equal to 1 does not correspond to any 
embedding of graphs on an oriented surface. Actually the definition of x can 
be generalized to (embeddings of graphs on) nonorientable closed surfaces £ by 
introducing the so-called "crosscap number" c, namely the number of Mobius 
bands one must attach to the sphere S 2 to obtain the desired surface. Such a 
"nonorientable genus" is related to the Euler number of S c by x(£ c ) — (2 — c) 
and x = 1 corresponds to the (minimum) crosscap number c = 1 characterizing 
the real projective plane P (MP 2 in the standard mathematical notation). 

The 9j(I) diagram was explicitly shown to be embeddable in P in 0, but 
this feature is actually well known because of its isomorphism with the complete 
bipartite graph K 3j3 (see e.g. IfTTTo R 

Taken for granted this result, the following Proposition deals with the existence 
of a well-defined insertion operation characterized by the the fact that it preserves 
the Euler number, namely the variation A n _i „ x of X m passing from 3(j — 1)(I) 
to 3nj(l) diagrams is zero. 

Proposition 4. 

There exists (up to symmetry) one kind of "twisted" insertion operation, denoted 
by Utx,, that generates the 3nj(I) from the 3(n — l)j(I) for any n > 3 according 
to the recurrence formula (0)- Moreover 3 M is topology-preserving with respect 
to embeddings in P. 

The proof can be carried out by resorting to the graphical method as already 
done in section 3.1 for type II diagrams. In the figure below the generation of the 
3nj(I) from the 3(n — l)j(I) according to the algebraic recurrence formula £[]) is 
shown. The summation over x implies the cancellation of both edges labeled by 
x and nodes which lose an incident edge; the last step would consist in joining 
back edges with the same labels to get the 3nj(T) diagram represented at the be- 
ginning of this section. Note that, up to relabeling, the same construction could 
be performed on anyone of the (n — 1) rays of the 3(n — l)j(I) and the resulting 
coefficients have the same values owing to their symmetry properties 0. 

3 At the end of section 2.1 the embedding of 9j(T) into the orientable, genus-1 torus T was 
considered. The fact that ^3,3 has the same nonorientable and orientable genus is exceptional: 
most graphs have different values for g and c. 
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The diagrammatic [algebraic] form of the twisted insertion operator J M to be used 
recursively for the generation of any 3nj(X) diagram [coefficient] is summarized 
in the following table (where the factor (2x + 1) and a phase have been dropped) 




n-l 



» 



31 



"II— 1 



kfl— 1 
/ A- 



ln—1 ^"n jn 



The topological invariance of such an operation is easily proved by observing that 
the variation of the Euler number owing to one insertion is given by 



A^xPm] = Ay -AE + AF = 2-3 + 1 = 



(14) 



where AV is the number of new vertices, AE the number of new edges and AF 
the number of new faces. (Recall that any type of 3nj diagram has 2n vertices 
and 3n edges and note that the (n + 1) "faces" for type I can by quickly selected 
by looking at the n cycles that include each of the n rays and adding +1 for the 
external region.) 

Now the proof of Proposition 3 can be completed by induction. Given that 
X [9j] = 1 and assuming x [3(n — = 1, by ([141) we get 



X [3nj] = x[3(n - l)j] + A n _i )7l xP> 



(15) 



for any n > 3. 
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A comment is in order about the possibility of generating recursively the other 
types of 3nj diagrams arising for n > 4 [|2l[H- We argue that, on applying suitable 
sequences of N insertion operations 3u (Prop. 2) and 3^ (Prop. 4) starting from 
the simplest graph -namely the complete quadrilateral associated with the 6j- all 
types of 3(2 + N)j diagrams can be generated. Actually the topological invariance 
proved in the two propositions holds true only when 3 n is applied to a type II 
diagram (and similarly 3^ to a type I diagram) under the circumstances described 
above. However, different positions of the insertions and/or actions on the other 
type of diagrams would provide a variety of configurations with possible changes 
of the underlying topology. Work is in progress in this direction. 

In view of applications to be addressed in the following we conclude this sec- 
tion by noticing that every 3nj(I) diagram has a Hamiltonian circuit of length 2n, 
as the samples of 12j, 15j, 18j(I) shown below make it manifest. The Hamilto- 
nian cycle would be always consistently labeled by the ordered set . . . , j n , 
hi, k 2 , ■ ■ ■ , k n and the rays by l u l 2 , . . . , in- 




Other non-Hamiltonian cycles appear as well in these diagrams. There are both 
cycles of length equal to n + 1 (including one ray and n edges on the perimeter) as 
well as cycles of length equal to the girth number, which is 4 for each n as already 
happened for 3nj(H) diagrams (now these 4-cycles are associated with "twisted" 
configurations of four edges labeled cyclically by hi, ji, k). 



4 Asymptotic disentangling of 3nj diagrams 

In two previous papers Il20l I2TI the mathematical apparatus of angular momenta 
recoupling theory 0] (topic 12) has been shown to provide a unifying background 
structure for a number of phenomena and applications ranging from atomic, molec- 
ular and nuclear physics to quantum computing and gravity (cfr. also Il2~2~l0 . The 
term spin networks, originally introduced by Penrose G3ll in connection with a 
model combinatorial spacetime, was used as a synonym of 3nj diagrams, namely 
connected and closed regular graphs whose 2n nodes are associated with triads 
of angular momentum variables, each satisfying triangle inequality. By semiclas- 
sical spin networks we mean, as in [120112711 . those families of 3nj diagrams that 



19 



arise when all or just a few of the 3n spin variables become large (namely the 
angular momentum labels attached to them are ^> 1 in h units). The issue of 
semiclassical analysis of angular momentum quantum transition amplitudes dates 
back to the beginning of quantum mechanics and we do not insist here on the wide 
range of applications of such techniques. Rather, we are going to reach a general 
assessment starting from what has been done for the 9j symbol in IfTOll . where the 
case in which six edges become large was shown to provide the disentangling of 
the underlying network, in the sense explained in section 4.1 below. (Note that 
a more formal approach based on geometric quantization techniques for both 3j 
(Wigner) coefficient ll24l and 6j symbol lfT3l can be carried out as well.) 

Whatever kind of asymptotics we want to address, the first object to be an- 
alyzed is the 3-valent node itself, namely the basic building block of any 3nj 
diagram. Given that a node is naturally associated with a 3j (Wigner) coefficient, 
and forgetting its dependence on magnetic quantum numbers, the possible choices 
of small/large angular momentum values are summarized in the following picture, 
where solid lines and small letters stand for truly "quantum" labels while capital 
letters and double dotted lines represent "semiclassical" large values. 



(quantum 3j) 




.A B. 



A B. 



C 



(0, 3) case (2, 1) case (3, 0) case 

Here we recover graphically the well known admissible asymptotics for 3j coeffi- 
cients, that is "if one spin value is large then at least one of the other incident spins 
must be large as well". Moreover, if the right-hand case (3, 0) occurs at anyone 
of the 2n node of a 3nj diagram, the resulting configuration turns out to be related 
to an amplitude whose distribution represents (a suitable number of) decoupled 
harmonic oscillators. Thus the "quantum entanglement of states", measured pre- 
cisely by the transition amplitudes associated with 3nj coefficients, is lost when 
approaching the classical limit, namely when the underlying quantum system un- 
dergoes a complete decoherence. 

The prototype of this phenomenon was addressed in the famous paper by Pon- 
zano and Regge [7] on "Semiclassical limit of Racah coefficients" dealing with 
the simplest spin network, the 6j symbol itself. In what follows we are going 
to explore and classify asymptotic regimes in which different kinds of "asymp- 
totic disentangling" arise in configurations that include (2, l)-nodes. (Recall that 
the well-known formulas for the (3, 3) and (4, 2) 6j cases, corresponding to the 
cycles shown in section 2.1, are listed as Eqns. 2 and 3 in ifTOlO . 
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As a case study we review first the main results on the 9j(I) proved in IfTOll by 
supporting them with new numerical simulations. 

4.1 The 9j case 

The focus of this section is on the 9j(I) symbol, well known in atomic spec- 
troscopy for its role as the matrix element of the transformation between LS and 
jj coupling schemes, but exhibiting features prototypical of more complex spin 
networks. In quantum chemistry, it appears for example in matrix elements when 
Sturmian or momentum space orbitals are employed [|25ll . On the basis of the 
remarks at the end of section 2.2, it is worth noting that the physical situations 
where some of the angular momenta become large, correspond mathematically to 
the case of discrete functions becoming continuous in the limit, as we are going to 
check in what follows. Note however that we are not going to address the analog 
of Ponzano-Regge asymptotics (9 entries of the 9j large, or (9, 0) expansion) and 
refer the reader to [5] Sec. 10.7 for details (the same reference is the basic one for 
notations and algebraic expressions used in this section). 

All cases where some entries are large and some small can be reduced to the 
(4, 5) or "large corner" case and to the (6, 3) or "small diagonal" case by permu- 
tation symmetries. The associated Yutsis graphs with dotted edges representing 
large entries are 




(4, 5) case (6, 3) case 

Recall that the basic expression of the 9j as a sum over x of the product of three 
6j was already given in (HJ). By slightly changing spin labels to comply with 
notations in ifTOl. it is rewritten as 




where max{|a — i\,\b — f\,\h — d\} < x < min{a + i,b + f, h + d}. 
The (6, 3) asymptotic case is derived from this defining equation by suitably ap- 
plying the well known (3, 3) expansions for the Qj [|71[8l and the final result reads 
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(recall that capital letters stand for large entries) 




B+H-e 

- d c-B,G-D (cos#) d* F _ CH _ G (cos9), (17) 



ri r 2 

where there appear two Jacobi polynomials (Wigner d-matrices of a cosine), with 
principal quantum numbers a, i (two of the "small" entries) and argument given 
by 

(2e + l) 2 - (B + F + l) 2 - (D + H + l) 2 



cos 9 



2(B + F + 1)(D + H + 1) 



Note that the role of magnetic quantum numbers in d J MM , is played here by dif- 
ferences of "large" spin variables which vary consistently between —j and +j in 
integer steps. Finally, the multiplicative factors T read 

T t = {ll(B + C + f) + l}ll(D + G + H) + l]} 1/2 ; 

T 2 = [(B + F + 1)(D + H + 1)] 1/2 . 

The (4, 5) case was addressed years ago by Watson [26] and his results are sum- 
marized in 




\b+C-d-G 



b e h \ ( d e f 



2X+1 \C-A e G-I)\G-A e C-I 

(18) 

where there appear two Wigner 3j symbols and X = \(A + C + G + /), e = 
A-C-G + I. 

Both formulas can then be taken as an illustration of "disentangling" of the as- 
sociated 9j networks because in the above semiclassical limits no summation ("su- 
perposition") appears anymore. Both d-matrices and Clebsch-Gordan (Wigner) 
coefficients represent well defined classes of orthogonal polynomials belonging to 
the Askey hierarchy as illustrated for instance in [1~3|| in connection with the quan- 
tum theory of angular momentum. As is well known, the issue of separability of 
variables in quantum mechanics is associated, by group theoretical arguments, to 
continuous or discrete symmetries, and allows the expression of wavefunctions of 
a composite system as a product of wavefunctions of subsystems. On the other 
hand, nonseparability requires introduction of coupling by expanding on a lin- 
ear combination of a set of such basis wavefunctions. Interpreting the original 
expression of the 9j given in © (or in (fT6l) ) as an expansion of a global "spin net- 
work" wavefunction, where the entries are discrete variables, the asymptotics (TT71) 
and (fT8~l) actually lead to wavefunctions proper of separated systems (note that the 
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surviving angular momentum functions depend on both truly quantum discrete 
variables and on semiclassical, continuous ones). Indeed, from the analysis car- 
ried out in section 2.1, we see that the variables are explicitly designated by c, d 
and h in the discrete case of a coupled set of difference equations, and by x, y, z in 
the continuous limit of a coupled set of partial differential equations. This appears 
to have been anticipated in an early remark by Neville lfT2~ll . 

From the graphical viewpoint, what is left after "deletion" of the dotted edges 
-associated to real-valued variables that tend to infinity- is in both cases a graph 
with no more closed cycles and thus "trivial". They look like isolated edges (case 
(6, 3)) or tree configurations (case (4, 5)), which are both topologically "con- 
tractible". The "entangled" nature of higher 3nj cases can be discussed along 
similar lines and according to remark (ii) of section 2.2. In the next section we 
are going to deal systematically with their disentangling. Note finally that such 
disentangling phenomena might be discussed and interpreted from many other 
viewpoints since spin network graphs (3nj diagrams) are encountered in quan- 
tum gravity models, discrete mathematics, and quantum computing ( 11201 l2il 1221 
and reference therein). We plan to investigate new applications of our asymptotic 
techniques in these contexts. 

Further important insight on disentangling may come from computational stud- 
ies. In ifTOl Wigner d-functions, Clebsch-Gordan coefficients and other special 
functions were calculated by directly summing their defining series using multi- 
precision arithmetic (MPFUN90) [|27ll . The multiple precision arithmetic allows 
convenient calculation of hypergeometric functions, of small and large argument 
by their series definition, and without the necessity of using recurence relations, 
integral and rational representations, or asymptotic approximations. On the basis 
of defininig relations and formulas for the functions given in various chapters of 
the basic problem encountered in these summations is that some of the terms 
in the sum may become very large (~ 10 100 or larger). However the sum must 
be a number in a small finite range, say [—1,1], and normal computer precision 
(8 to 32 decimal places) will not generally produce accurate summation. On the 
other hand, multiple precision arithmetic packages that allow calculations with 
thousands (or millions!) of digits sum the hypergeometric series easily for a spin 
variable j of magnitude much larger than one thousand. The asymptotic formulas 
(fT71) and (PT8l have been already tested in IfTOll by calculating 9j symbols and their 
approximations for huge randomly selected sets of j values but we conclude this 
section by presenting new numerical experiments, illustrated in the following and 
displayed in the figures below. 

The large entries on the coefficients (either 6 or 4 according to either (fT71) or 
(fT8l) ) have values ranging from a given abscissa, J to J + j m , where J is chosen 
between 20 to 100. The small j (integer and half integer values) are randomly 
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selected from the range [0, j m ]. The calculations use 1000 nonzero 9j values for 
each J and integer and half integer small j for j m = 2, 9/2, 19/2, 29/2. We cal- 
culate the rms deviations, i.e. the root mean square deviations of the approximate 
9j from the exact values, and the rms magnitudes of the 9j symbols. We take 
the ratio of these quantities to evaluate the fractional error in the disentangled 9j 
symbol. Then the numerical experiments show that the fractional error decreases 
as J -1 for both Equations (fT71) and (PT8l) ) while the rms magnitudes decrease as 
J- 2 for Eq.dn]) and J- 1 for Eq.(Q]§). 

Experiments with different values for j m show that the fractional errors for 
both equations scale as j^J J. For the six large angular momenta formula, the 
fractional error is approximately 0.3j,^/J. For the four large angular momenta 
formula, the fractional error is approximately 0.2j^/ J . The Watson equation (PT8l) 
produces a significant number of zero values (15 to 20 %) for small but nonzero 
exact 9j symbols. These zeros result from combinations of angular momenta 
that yield zero values for the 3j symbols. Preliminary numerical checks indicate 
a correlation with the previously mentioned phenomenon of the appearance of 
nodes in separable conditions. Remarkably the rms magnitudes of the 9j symbols 
that correspond to zero approximate values scale as J~ 2 and have approximately 
the same values as the rms deviation values. Hence the quality of the overall 
approximation not degraded even if the cases with zero approximate values are 
included. 
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Figure 1: 
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Figure 3: 
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4.2 The general case 

Recall from section 3 that both 3nj(T) and 3nj(II) spin networks possess a Hamil- 
tonian circuit of length 2n. In this section we are going to show how disentan- 
gling -induced by letting the 2n spin variables on this closed cycle to become 
large- works for any 3nj diagram of type I by applying recursively an asymptotic 
version of the insertion operation 3^ (a similar procedure might be set up for 
3nj(U)). 

In order to make the notation used in section 4.1 compatible with the general 
case, let us rewrite the 9j(I) accordingly in terms of the labels {j, k, 1} (see © and 
the diagram of the 3nj(T) at the beginning of section 3.2). Consider at the same 
time the limit in which the Hamiltonian cycle (the external circle of the cartwheel 
diagram with 3 internal edges) is labeled by large spin values, denoted by capital 
letters 



3i J2 J3 . , , , , 
H h f 3 > > 



k\ ki 



J\ J2 J3 

h k k > (19) 

K x K 2 K 3 



To such a (6, 3) expansion of the 9j (6 large, 3 small entries) we can apply the gen- 
eral twisted insertion operation defined in Prop. 4 of section 4.2 to get the 12j(I). 
In the present case, without imposing any restriction on the entries, it would read 





k 3 x 1 
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X 1 


ft 


£4 fc 4 J 


U 


k 


3 a] 



(20) 



However, since such a combination has to be "coupled" with the (6, 3) 9j (in 
which an x takes formally the place of the "small" l 3 in ( fl9l ) we see that the entries 
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ji, k 3 , k±, jz are necessarily "large" because they lie on the preexisting Hamilto- 
nian great circle, while Z 3 stays small. Then we are left with the residual problem 
of determining the order of magnitude of the new entries j 4 , fc 4 , / 4 (recall that x is 
small). This is achieved by checking the triads of the two 6j symbols in (1201) to 
see which possibilities we have for these three labels in view of the admissibile 
configurations for the 3j coefficient discussed in the introduction to this section. 

. The triads of { £ £ * } are 

( Ji K 3 x) (2 large, 1 small) and as such it is admissible; 
(I3 I4 x) (2 small) and thus Z 4 must be small; 

( Ji U fc 4 ) (1 large, 1 small) and then fc 4 must be large and renamed K4, 
(l 3 K 3 K4) (2 large, 1 small) and then admissible. 

. The triads of { § J * } are 

(Ki J 3 x) (2 large, 1 small) and as such it is admissible; 
(l 3 I4 x) is the same triad considered before; 

(Ki (1 large, 1 small) and then j 4 must be large and renamed J 4 ; 
(h J3 Ja) (2 large, 1 small) and then admissible. 

Since the above remarks clearly hold true for any insertion operation which 
acts on a preexisting [2(n — 1) large, (n — 1) small] 3(n — l)j(I) diagram, we have 
established the following 

Proposition 5. 

There exists (up to symmetry) one kind of twisted "asymptotic" insertion opera- 
tion, denoted by 3^ , that generates a unique Hamiltonian disentangled (2n, n) 
3nj(I) diagram from the similar [2{n — 1), (n — 1)] 3(n — l)j(I) for any n > 3. 

Its general (graphical and algebraic) expression reads 



, : Ji 
•J 



J, 



11- 



J\ Kn-i x I J J\ Jn-i x 



• hi \ «!, 

where edges with large labels are drawn as dotted lines. 

Thus we can easily infer that, on the basis of the (6, 3) expansion of the 9j 
given in (TT71) . the analytical expression of any Hamiltonian asymptotic expansion 
of a (2n,n) 3nj(T) diagram contain the product (no sum) of (n — 1) Wigner d- 
functions whose principal quantum numbers may be chosen among the n spin 
variables Z 2 , • • • , l n } (a similar result can be established for 3nj(ll) diagrams). 
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5 Conclusions and outlook 



We have developed in this paper a recursive procedure for the generation of 3nj 
from 3(n — \)j diagrams of types I and II as well as an original approach to the 
classification of their asymptotic disentangling. 

There are a number of open questions that deserve further investigations. 
The recursive generation of all other types of 3nj diagrams would be a major 
achievement, as pointed out already in [Q]| (topic 12). We argue that it would be 
possible (at least in principle) to complete the enumeration of diagrams of type III, 
IV, V,... for any n > 4 starting from either type I or II and applying suitably chosen 
sequences of the two kind of insertion operators defined in section 3. However it 
should be clear that the enumeration will blow up as n increases and the hardest 
task will be to rule out isomorphic configurations (see footnote 2 and IIT810 . 

For what concerns semiclassical limits and asymptotic disentangling, we are 
currently performing numerical experiments on the five types of 15j coefficients 
to test the relative magnitudes and associated probabilities (cfr. If2~8l for an early 
attempt of comparison between type I and II). We have restricted most of the 
treatment in section 4 to Hamiltonian circuits, but it can be conjectured that the 
described features of disentangling will occur for other circuits in general spin 
networks. It would be also interesting to analyze the different asymptotic regimes 
in connection with the symmetry properties of the different types of coefficients. 
An improvement of multi-variable recursion relations and related partial differ- 
ential equations (see the 9j case discussed in section 2.2) would greatly help in 
providing further insights in the classification of the associated orthogonal poly- 
nomials of hypergeometric type. 

The results obtained here, as well as possible new improvements along the 
lines outlined above, represent a prerequisite for extending to 3nj the analysis 
based on sophisticated (quantum and semiclassical) geometric and analytical tech- 
niques introduced in [|TTll24l . see also 11291 . 
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